Post-Newtonian Approximation in Maxwell-Like Form 
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The equations of tlie linearized first post-Newtonian approximation to general relativity are of- 
ten written in "gravitoelectromagnetic" Maxwell-like form, since that facilitates physical intuition. 
Damour, Soffel and Xu (DSX) (as a side issue in their complex but elegant papers on relativistic 
celestial mechanics) have expressed the first post-Newtonian approximation, including all nonlin- 
earities, in Maxwell-like form. This paper summarizes that DSX Maxwell-like formalism (which is 
not easily extracted from their celestial mechanics papers), and then extends it to include the post- 
Newtonian (Landau-Lifshitz-based) gravitational momentum density, momentum flux (i.e. gravita- 
tional stress tensor) and law of momentum conservation in Maxwell-like form. The authors and their 
colleagues have found these Maxwell-like momentum tools useful for developing physical intuition 
into numerical-relativity simulations of compact binaries with spin. 



I. INTRODUCTION 

In 1961, Robert L. Forward [H (building on earlier 
work of Einstein 0, (3 and especially Thirring 0, 
wrote the linearized, slow-motion approximation to gen- 
eral relativity in a form that closely resembles Maxwell's 
equations; and he displayed this formulation's great in- 
tuitive and computational power. In the half century 
since then, this Maxwell-like formulation and variants 
of it have been widely explored and used; see, e.g., 
d, 0, S i, M, HH, m, [13 and references therein. 

In 1965-69 S. Chandrasekhar [H, ^ formulated the 
first post-Newtonian (weak-gravity, slow-motion) approx- 
imation to general relativity in a manner that has been 
widely used for astrophysical calculations during the 
subsequent 40 years. When linearized, this first post- 
Newtonian (IPN) approximation can be (and often is) 
recast in Maxwell-like form. 

In 1991, T. Damour, M. Soffel and C. Xu (DSX; 0) 
extended this Maxwell-like IPN formalism to include all 
IPN nonlinearities (see also Sec. 13 of Jantzen, Carini 
and Bini [13). DSX did so as a tool in developing a 
general formalism for the celestial mechanics of bodies 
that have arbitrary internal structures and correspond- 
ingly have external gravitational fields characterized by 
two infinite sets of multipole moments. (For a generaliza- 
tion to scalar-tensor theories, see [13 ■) In 2004, Racine 
and Flanagan [l9!| generalized DSX to a system of com- 
pact bodies (e.g. black holes) that have arbitrarily strong 
internal gravity. 

During the past 18 months, we and our colleagues have 
been exploring the flow of gravitational field momentum 
in numerical-relativity simulations of compact, spinning 
binaries j20l. \2]\. In their inspiral phase, these binaries' 
motions and precessions can be described by the IPN ap- 
proximation,^ and we have gained much insight into their 



dynamics by using the IPN DSX Maxwell-like formalism, 
extended to include Maxwell- like momentum density, mo- 
mentum flux, and momentum conservation.^ 

In this paper, we present that extension of DSX,'^ 
though with two specializations: (i) we fix our coordi- 
nates (gauge) to be fully harmonic instead of maintain- 
ing the partial gauge invariance of DSX, and (ii) we dis- 
card all multipole moments of the binaries' bodies except 
their masses and their spin angular momenta, because for 
black holes and neutron stars, the influences of all other 
moments are numerically much smaller than 1.5PN or- 
der. 

The DSX celestial-mechanics papers [H, [H [H [111 
are so long and complex that it is not easy to extract 
from them the bare essentials of the DSX Maxwell-like 
IPN formalism. For the benefit of researchers who want 
those bare essentials and want to see how they are related 
to more conventional approaches to IPN theory, we sum- 
marize them before presenting our momentum extension, 
and we do so for a general stress-energy tensor, for a per- 
fect fluid, and for a system of compact bodies described 
by their masses and spins. 

This paper is organized as follows: In Sec. |TT] we 
summarize the basic DSX equations for IPN theory in 
Maxwell-like form. In Sec. IIIII we specialize the DSX 
formalism to a perfect fluid and make contact with the 



^ For black-hole and neutron-star binaries, the influences of spin 
that interest us are formally IPN, but because of the bodies' 
compactness (size of order Schwarzschild radius), they are nu- 



merically 1.5PN. 
^ A referee has pointed out to us that some papers in the rich liter- 
ature on the Maxwell-like formulation of linearized IPN theory, 
e.g. [10| . argue that the Maxwell analogy is physically useful only 
for stationary phenomena. Our spinning-binary application [201 
of the momentum-generalized DSX formalism is a counterexam- 
ple. 

' When we carried out our analysis and wrote it up in the 
original version of this paper, we were unaware of the 
Maxwell-like formalism in DSX 'l(jfi : see our preprint at 
'http://xxx.laiil.gov/abs/0808.2510vl When we learned of 
DSX from Luc Blanchet, we used it to improve our Maxwell-like 
treatment of gravitational momentum (by replacing our defini- 
tion for the gravitoelectric field by that of DSX) and we rewrote 
this paper to highlight the connection to DSX. 
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conventional IPN notation. In Sec. IIVI we extend DSX 
by deriving the (Landau-Lifshitz-based) density and flux 
of gravitational momentum in terms of the DSX gravito- 
electric and gravitomagnetic fields and by writing down 
the law of momentum conservation in terms of them. (It 
is this that we have found so useful for gaining intuition 
into numerical-relativity simulations of inspiraling, spin- 
ning binaries 20,21].) In Sec. |V] we briefly discuss energy 
conservation. In Sec. IVIl relying on Racine and Flana- 
gan [l^ . we specialize to the vacuum in the near zone 
of a system made from compact bodies with arbitrarily 
strong internal gravity. Finally, in Sec. IVIII we summa- 
rize the DSX formalism and our extension of it both for a 
self-gravitating fluid and for a system of compact bodies. 

Throughout this paper, we set G = c = 1, Greek letters 
run from to 3 (spacetime) and Latin from 1 to 3 (space) , 
and we use the notation of field theory in flat space in a 
34-1 split, so spatial indices are placed up or down equiv- 
alently and repeated spatial indices are summed whether 
up or down or mixed. We use bold-face italic characters 
to represent spatial vectors, i.e. w is the bold-face version 

of Wj. 



II. THE DSX MAXWELL-LIKE FORMULATION 
OF IPN THEORY 

Damour, Soffel and Xu (DSX [H]) express the IPN 
metric in terms of two gravitational potentials, a scalar 
w and a vector Wj : 

500 = -e-^"" = -l + 2w-2w^+0{Ufj) , 

go^ = -4u;, + 0([/^/') , (2.1) 

g,^ = <5^,e2» = 5„(1 + 2w) + 0{Uj,) 

[DSX Eqs. (3.3)]. The Newtonian hmit of w is Un ^ 
(Newtonian gravitational potential), and Wi is of order 



U 



3/2 
N • 



w = Un + 0{Uf,) , = 0{ulj^) . (2.2) 

The harmonic gauge condition implies that 

w.t + Wj,j = (2-3) 

[DSX Eq. (3.17a)]; here and throughout commas denote 
partial derivatives. Using this gauge condition (which 
DSX do not impose), the IPN Einstein fleld equations 
take the following remarkably simple form: 

V^w -w = -4tt{T"° + T") + 0{U%/C^) , (2.4a) 
V^w, = -AttT^' + 0{Ull^ / C^) (2.4b) 

[DSX Eqs. (3.11)]. Here T"'' is the stress-energy ten- 
sor of the source (which we specialize below to a perfect 
fluid), is the flat-space Laplacian (i.e. V^w = Wjj), 
repeated indices are summed, dots denote time deriva- 
tives (i.e. w = vj^tt), and C is the lengthscale on which w 
varies. 



Following DSX, we introduce the IPN gravitoelectric 
field g (denoted e or by DSX , depending on the con- 
text) and gravitomagnetic field H (denoted 6 or B by 
DSX): 

g = Vw + Aw + OiUls/C), (2.5a) 
H = -4V XW + 0{U^^/C) . (2.5b) 

[DSX Eqs. (3.21)]. 

The Einstein equations (|2.4p and these definitions im- 
ply the following IPN Maxwell-like equations for g and 
H: 



Vg 

V X gr 

V H 
V xH 



-Att{T°" + T^^)~3w■ 
H + 0{gUl/C) , 
Q + 0{HUn/C) , 



0(gC/^/£X?-6a) 
(2.6b) 
(2.6c) 

167rr°Vj +Ag + 0{HUnIC) (2.6d) 



[DSX Eqs. (3.22)]. Here is the unit vector in the i 
direction. 

In terms of g and H , the geodesic equation for a parti- 
cle with ordinary velocity v — dx/dt takes the following 
form [Eq. (7.17) of DSX, though in a less transparently 
"Lorentz-force" -like form there] : 



^[(1 
dt '■^ 



3[/jv 



v] = {1-Un + lv'^)g + v X H 
+0{gU^^) . (2.7) 



Note that the spatial part of the particle's 4- 
momentum is mu — m(l + Un + ^v^)v at IPN order. 
This accounts for the coefficient 1 -|- Un + on the 
left-hand side of Eq. (|2.7p . The remaining factor 2Um 
is related to the difference between physical lengths and 
times, and proper lengths and times. In the linearized, 
very-low-velocity approximation, this geodesic equation 
takes the "Lorentz-force" form dv/dt = g + v x H, first 
deduced (so far as we knowHn 1918 by Thirring [3], mo- 
tivated by Einstein's 1913 2] insights about similarities 
between electromagnetic theory and his not-yet-perfected 
general relativity theory. 

The IPN deviations of the geodesic equation (|2.7|) from 
the usual Lorentz-force form might make one wonder 
about the efficacy of the DSX definition of g. That ef- 
ficacy will show up most strongly when we explore the 
gravitational momentum density in Sec. IIVI below. 



III. SPECIALIZATION TO A PERFECT FLUID 

We now depart from DSX by specializing our source to 
a perfect fiuid and making contact with a set of IPN grav- 
itational potentials that are widely used. We pay special 
attention to connections with a paper by Pati and Will 
[13] because that paper will be our foundation, in Sec. 
IVI for computing the density and flux of gravitational 
field momentum. 
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We describe our perfect fluid in the following stan- 
dard notation: po — (density of rest mass), 11 = (internal 
energy per unit rest mass, i.e. specific internal energy), 
P = (pressure), all as measured in the fluid's local rest 
frame; Vj = dx^ /dt — (fluid's coordinate velocity). 

Following Blanchet and Damour [28| . and subsequently 
Pati and Will (Eqs. (4.13), (4.3) of [H), we introduce a 
post-Newtonian variant U of the Newtonian potential, 
which is sourced by + T": V^C/ = -47r(T °" + T "). 
Accurate to IPN order, the source is [see Eq. (|4.4dp ] 

7^00 ^ ^ + Il + 2v^ + 2Un + iP/po) , (3.1) 

where Un is the Newtonian limit of U 



(3.2a) 



Correspondingly, U can be written as 



U 



Po{l + Il + 2v^ + 2Un + iP/Po) 



\x — X' 



d-'x' . (3.2b) 



Here and below the fluid variables and gravitational po- 
tentials in the integrand are functions of {x',t) as in Eq. 
p.2ap . In Eq. (|3.2ap for Un, po can be replaced by any 
quantity that agrees with po in the Newtonian limit, e.g. 
by the post-Newtonian "conserved mass density" of 
Eq. (j3.6bp below. We also introduce Chandrasekhar's 
Post-Newtonian scalar gravitational potential x (Eq. (44) 
of [III), which is sourced by 2Un, V^x = -2Un or equiv- 
alently 



X 



J Po\x - x'\Sx' . 



(3.3) 



Pati and Will use the notation —X for x (Eqs. (4.14), 
(4.12a) and (4.3) of [13). 

It is straightforward to show that the IPN solution to 
the wave equation (|2.4ap for the DSX scalar potential w 
is 



w 



U 



1 



(3.4) 



and the IPN solution to the Laplace equation (|2.4bp for 
the DSX vector potential Wj is 



-^^d^x' 



(3.5) 



The fluid's evolution is governed by rest-mass conserva- 
tion, momentum conservation, and energy conservation. 

The IPN version of rest-mass conservation takes the 
following form: 



p*,t + V • (p^v) = , 



where 



PoU 



Po{l 



1 2 
1^ 



3C/) 



(3.6a) 



(3.6b) 



(Eqs. (117) and (118) of Chandrasekhar fli|). Here u° is 
the time component of the fluid's 4-velocity and g is the 
determinant of the covariant components of the metric. 

We shall discuss momentum conservation and energy 
conservation in the next two sections. 

We note in passing that Chandrasekhar and many 
other researchers write their IPN spacetime metric in a 
different gauge from our harmonic one. The two gauges 
are related by a change of time coordinate 



tc — ta — ^x 



(3.7a) 



and correspondingly the metric components in the two 
gauges are related by 



500 ^ 9oo 



(3.7b) 



Here C refers to the Chandrasekhar gauge and H to our 
harmonic gauge. DSX write their equations in forms that 
are invariant under the gauge change (j3.7p . 



IV. MOMENTUM DENSITY, FLUX, AND 
CONSERVATION 

We now turn to our extension of the DSX formalism to 
include a Maxwell-like formulation of gravitational mo- 
mentum density, momentum flux, and momentum con- 
servation. Following Chandrasekhar [l5[, Pati and Will 
[131 and others, we adopt the Landau-Lifshitz pseudoten- 
sor as our tool for formulating these concepts. 

From Pati and Will's 2PN harmonic-gauge Eqs. (2.6), 
(4.4b) and (4.4c) for the pseudotensor, one can deduce 
the following IPN expressions for the the gravitational 
momentum density and momentum flux (stress) in terms 
of the DSX gravitoelectric and gravitomagnetic fields g 
and H: 



(-.9)^LLei 



(4.1a) 
(4.1b) 



Each equation is accurate up to corrections of order Un 
times the smallest term on the right side (2PN correc- 
tions). 

For comparison, in flat spacetime the electromagnetic 
momentum density is j^E x B and the momentum flux 
is ^{E.E, - lS,,EkEk) + i^{B,Bj - \5,jBkBk). Aside 
from a sign in Eq. (^TTop and the two terms involving 
Un , the gravitational momentum flux and density j^. J[ ) 
are identical to the electromagnetic ones with E ^ g 
and B ^ H . Therefore, by analogy with the electro- 
magnetic case, there are gravitational tensions |grp/87r 
and |iifp/87r parallel to gravitoelectric and gravitomag- 
netic field lines, and gravitational pressures of this same 
magnitude orthogonal to the field lines. This makes the 
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gravitoelectric and gravitomagnetic fields g and H pow- 
erful tools for building up physical intuition about the 
distribution and flow of gravitational momentum. We 
use them for that in our studies of compact binaries [13] , 
relying heavily on Eqs. (|4.ip . 

Here are some hints for deducing Eqs. (|4.1[) from Pati 
and Will [IJ (henceforth PW): (i) Show that the last two 
terms in (2.6) of PW are of 2PN order for {a, (3} = {Qj} 
or {ij} and so can be ignored, whence lQ'K{—g)t'^^ — 
K°'^ . (ii) Show that our notation is related to that of 
PWby X = -X, [/ the same, w = U-\x = \{N + B)- 
\{N + Bf [for the last of these cf. PW (5.2), (5.4a,c)], 
and at Newtonian order Un = \N . (iii) In PW (4.4b,c) 
for A"^, keep only the Newtonian and IPN terms: the 
first curly bracket in (4.4b) and first and second curly 
brackets in (4.4c). Rearrange those terms so they involve 
only K, N+B and the Newtonian-order N, use the above 
translation of notation and use the definitions (|2.5p of g 
and H. Thereby, bring PW (4.4b,c) into the form (|4.ip . 

In the Landau-Lifshitz formalism, the local law of 4- 



momentum conservation T-"^. 



takes the form 



(4.2) 



(Eqs. (20.23a) and (20.19) of [29|, or (100.8) of [30|). 
Here (as usual), commas denote partial derivatives, and 
semicolons denote covariant derivatives. This is the con- 
servation law that we use in our studies of momentum 
flow in compact binaries [20l |. 

When dealing with material bodies (e.g. in DSX) 
rather than with the vacuum outside compact bodies, 
an alternative Maxwell-like version of momentum con- 
servation is useful. Specifically, using expressions ()4.ip 
and {-g) = 1 + 4Un [from Eq. (|27T|) with w = Un at 
leading order], and using the field equations (12. 6p for g 
and H, the conservation law (|4.2p can be rewritten in 
the following simple Lorentz-force-like form (Eq. (4.3) of 
Damour, Soffel and Xu's Paper II 



[(1 + 4c/w)r" 



[(1 



4c/jv)r^] 



(4.3) 



{T°° + T")g, + e,,kT°^Hk 



Here the Levi-Civita tensor produces a cross prod- 
uct of the momentum density with the gravitomagnetic 
field. For comparison, in fiat spacetime, the momen- 
tum conservation law for a charged medium interacting 
with electric and magnetic fields Ei and Bi has the form 
T"°_t -I- T'-^ J = peEi + eykJjBk, where pe is the charge 
density and Jj the charge fiux (current density). The 
right-hand side of Eq. ( the gravitational force den- 
sity) is identical to that in the electromagnetic case, with 
Pe {T°° T«), Jj T°^, E ^ g, and B H. 
Again, this makes g and H powerful foundations for 
gravitational intuition. 

For a perfect fluid, the components of the IPN stress- 
energy tensor, which appear in the momentum conser- 
vation law (I4.3p . are (Eqs. (20) of Chandrasekhar [13]) 



r 



00 



Poil + Ii + v^ + 2Un - 
Po{l + Ii + v^ + 2Un 
+P{1 - 2Un)5.,j , 
/9(1 + H + 2v^ + 2Un 



(4.4a) 
P/po)«, (,4.4b) 

- P/po)ViVj 

(4.4c) 

- 3P/po) (4.4d) 



V. ENERGY CONSERVATION 

For a perfect fluid, the exact (not just IPN) law of en- 
ergy conservation, when combined with mass conserva- 
tion and momentum conservation, reduces to the first law 
of thermodynamics dH/dt = —Pd{l/po)/dt\ so whenever 
one needs to invoke energy conservation, the first law is 
the simplest way to do so. For this reason, and because 
deriving the explicit form of IPN energy conservation 
[[-~g)[T^^ + *ll)],m = is a very complex and delicate 
task (cf. Sec. VI of [3l|), we shall not write it down here. 

However, we do write down the Newtonian law of en- 
ergy conservation in harmonic gauge, since we will oc- 
casionally need it in our future papers. Chandrasekhar 
calculated {-g){T°^' + t^^) in his Eqs. (48) and 

(57) are the time-time and time-space components, re- 
spectively. When one writes the expressions in terms of 
the "conserved rest-mass density" [Eq. (j3.6bp ] and in 
our Maxwell- like form, Newtonian conservation of energy 
states that 

p*(i + n + + 3Un) ~-^a -9 

P 1 

p*t)(l + HH h-w^ + 3C/jv) 

p 2 



3 r-r 1 H- 

Att At: 



= 



(5.1) 



While this equation is perfectly correct, it expresses 
Newtonian energy conservation in terms of the post- 
Newtonian gravitomagnetic field H . It is possible to 
rewrite the i?-dependent term using the relationship, 
V • [-l/(4^)(flr xH)] = V ■ [-AUnT^^Bj + (l/7r)[/jvg], 
which is accurate up to corrections of order g ■ H . This 
relationship can be found by applying Eq. (j2.6bp once 
and (j2.6dp twice, in combination with elementary vector- 
calculus identities. The statement of Newtonian energy 
conservation then depends only upon the Newtonian po- 
tential and its gradient and time derivative: 



1 7 

p^i + H + + 3Un) -^g g 



p 1 

/9*t>(l + HH V -v^ - Un) 

p 2 



3 • 1 

--r^^g + -u^g 

Aw TT 



= 



(5.2) 
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Notice that going from Eq. (|5.ip to Eq. (|5.2p involves 
adding a divergence-free piece to the energy flux, so it 
entails changing how the energy flux is localized — a 
change that strictly speaking takes the energy flux out of 
harmonic gauge. 

If the coefficients of the gravitational terms in Eq. (|5.2p 
look unfamiliar, it is because even at Newtonian order, 
the density and flux of gravitational energy are gauge- 
dependent. In some other gauge, they will be different; 
see Box 12.3 of [H. 



VI. GRAVITATIONAL POTENTIALS IN THE 
VACUUM OF A SYSTEM OF COMPACT, 
SPINNING BODIES 

For a system of compact, spinning bodies (neutron 
stars or black holes), the gravitational potentials Un, U, 
Wi and X in the vacuum between the bodies take the fol- 
lowing forms (in a slightly different harmonic gauge than 
the one used in Sec. Illll for fluids): 



Un 



U 



E 



Ma 

TA 



(6.1a) 




-2E 



2^ A 




TA^ 



TA 



U--X, 



Mav\ 



E 



TA 



TA' 



(6.1b) 
(6.1c) 

(6.1d) 
(6.1e) 



VII. CONCLUSION 

In our Maxwell-like formulation of the IPN approxima- 
tion to general relativity for fluid bodies, the evolution of 
the fluid and gravitational fields is governed by: (i) the 
law of momentum conservation (|4.3[) , (|4.4p (which can be 
thought of as evolving the fluid velocity Vj); (ii) the law 
of mass conservation p.6p (which can be thought of as 
evolving the mass density po)] (iii) the equation of state 
P{po) and first law of thermodynamics dll — ~Pd{l/ po) 
(which determine P and 11 once po is known); Eqs. p.2p . 
(|3.5p . (|3.3p for the gravitational potentials U, x, and wj; 
and Eqs. (|2.5p or (|2.6p for the gravitoelectric and gravit- 
omagnetic fields g, H. 

When specialized to a system of compact bodies, e.g. 
a binary made of black holes or neutron stars, the sys- 
tem is governed by: (i) IPN equations of motion and 
precession for the binary (not given in this paper; see, 
e.g., Eqs. (4.10), (4.11) and (4.14) of [1^); (n) momen- 
tum flow within the binary as described by the Landau- 
Lifshitz pseudotensor (|4.ip and its conservation law (|4.2p , 
in which the gravitoelectric and gravitomagnetic fields 
are expressed as sums over the bodies via Eqs. (|6.1I) and 
(|2.5p : and (iii) other tools developed by Landau and Lif- 
shitz (Sec. 100 of [30|). We are finding this formalism 
powerful in gaining insight into compact binaries [2d| . 
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Here the notation is that of Sec. IV of Thorne and Hartle 
[26'| : the sum is over the compact bodies labeled by capi- 
tal Latin letters A, B; Ma, S^j^ and v\ are the mass, spin 
angular momentum and coordinate velocity of body A; 
rA is the coordinate distance from the field point to the 
center of mass of body A; tab is the coordinate distance 
between the centers of mass of bodies A and B; n\ is the 
unit vector pointing from the center of mass of body A 
to the field point; and e^fe is the Levi-Civita tensor. 

Equations (|6.ip for the potentials can be deduced by 
comparing our IPN spacetime metric coefficients [Eqs. 
(lO) . (lOl) ] with those in Eqs. (2.4), (5.11) and (5.14) of 
Racine and Flanagan [l^ or in Eqs. (4.4) of Thorne and 
Hartle [li].^ 



* Racine and Flanagan specialize DSX to a system of compact 
bodies with a complete set of nonzero multipole moments. We 



neglect all moments except the bodies' masses and spins (see 
fifth paragraph of Sec. |T]|. Our notation is related to that of 
Racine and Flanagan by Upf = w = U — = — (3> -|- 

ip),Wi = — jCi. The Racine-Flanagan derivation of Eqs. 116. Il l 
avoids considering the internal structures of the bodies and it 
therefore is directly valid for black holes. The Thorne-Hartle 
derivation relies on the pioneering analysis of Einstein, Infeld 
and Hoffman [33l | which uses a fluid description of the bodies' 
interiors. Thorne and Hartle extend Eqs. 1 16. Il l to black holes by 
th e eq uivalent of Damour's "effacement" considerations, Sec. 6.4 

of [13. 
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